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Beschreibung 

Verfahren zur Stcucrung der Pulssequenz fur eine Kernspintomographieanlago unci Vorrichtung zur Durchfiihrung des 
Verfahrens. * 
S Zum iMeBablauf bei einer K-jmr-p!niomogr:aph5C2*. , u2^ : : ::d ; ;r. v.'e;-cn:Ii ch-jn Zeitvcrlauf von Gradientenstromen, 
Hochfrequenz-Scndcpulse und Samplingperioden far das gewonnene MR-Signai zu sieusrn. Die enlsprechenden Sleu- 
erdatensatze wcrden bei den heute ublichen Steuerungen im wesentlichen vor dem Sequenzsiari vorberechnet. Dabei 
greifi man auf abgespeicherte Zeitverlaufe, z. B. fur die Anstiegs- und Abfallflanken, von Gradientenstromen zuriick. 
Wahrend des Sequenzablaufs wird dann der MeBvorgang in Form einer Liste "abgespielt". Eine derartige Steuerung ist 
10 z. B. aus der US-Patentschrift 5,349,296 bekannt. * 

Da der gesamte Sequenzablauf bei Sequenzstart bereits festgelegt ist, kann man nach dem Sequenzstart nur noch 
schwer und sehr eingeschrankt auf den weiteren Sequenzablauf EinfluB nehmen. Im wesentlichen ist nur ein eingeplantes 
Anhalten und Warten auf ein extemes Triggersignal moglich. Es muB ein groBes Datenvolumen abgespeichert werden. 
Um dieses Datenvolumen einigermaBen in Grenzen zu halten, werden typischerweise Teile der Sequenztopologie in 
15 Form einer Schleifenstruktur abgebildet. Im Datensatz sind also z. B. Schieifenanweisungen enthalten. Darnit wird die 
Beschreibungs- "Sprache" fur den Datensatz komplex und unflexibel. 

Aus der US-Patentschrift 5,144,242 ist eine Steuereinrichtung fur MR-Gerate bekannt, bei der ein Speicher fur Steu- 
erbefehle kontinuierlich aus einem Massenspeicher nachgeladen wird, wahrend die Sequenz abgearbeitet wird. Darnit 
wird der Speicher fur Steuerbefehle entlastet. Eine ahnliche Steuereinrichtung ist auch aus der US-Patentschrift 
20 5,606,258 bekannt. 

In der US-Patentschrift 5,512,825 ist ein Verfahren beschrieben, zur Minimierung von Totzeiten fur die Gradienten 
Anfangs- und Endwerte sowie ein Moment fur den Gradientenverlauf anzugeben. 

Bei der Definition von Gradienten ist zu unterscheiden zwischen den sogenannten "logischen" Gradienten und den 
"physikalischen" Gradienten. Diese Unterscheidung wird dadurch notwendig, daB mit Kemspintomographiegeraten be- 

25 liebig schrage Schichten abgebildet werden konnen. Diese schragen Schichten werden durch Vorgabe von logischen Gra- 
dienten definiert, die in einem kartesischen Koordinatensystem entsprechend schrag liegen. Jede Gradientenspule kann 
naturlich nur einen Gradienten in einer festgelegten Richtung, narnlich nur in einer Achse eines physikalischen kartesi- 
schen Koordinatensy stems liefern. Diese Gradienten werden als physikahsche Gradienten bezeichnet. Die schragen lo- 
gischen Gradienten werden daher durch Uberlagerung von physikabschen Gradienten erzeugt, mathematisch gesehen 

30 durch eine Vektoraddition. Darnit entsteht jedoch folgendes Problem: Bei herkommlichen Anlagen werden die logischen 
Gradienten typischerweise in Form von Trapez-Pulsen definiert. Bei einer Uberlagerung mehrerer logischer Gradienten 
sind die physikalischen Gradienten dann typischerweise nicht mehr trapezformig, sondern weisen einen komplexen 
Strom-Zeit-Verlauf auf, der vom Gradientenverstarker beherrscht werden muB. Der maximale Anstieg der Gradienten- 
flanken der physikalischen ist bei jedem Gradientensystem durch die sogenannte Slewrate begrenzt Femer ist naturlich 

35 auch die Gradienten- Amplitude begrenzt. Bei der Definition der logischen Gradienten wird nicht vorherbesummt, wel- 
che Anstiegszeiten und welche Amplituden sich bei den physikalischen Gradienten ergeben, da diese erst bei Sequenz- 
ablauf aus den logischen Gradienten berechnet werden. Andererseits muB aber bei der Definition von vorbereiteten Da- 
tensatzen fiir die logischen Gradienten die Begrenzung der Anstiegszeit und der Amplituden der physikalischen Gradien- 
ten berucksichtigt werden. Es muB z. B. immer vom ungtinstigsten Fall ausgegangen werden, daB alle Gradientenan- 

40 stiege auf einen physikalischen Gradienten trefFen. Darnit ist nur eine schlechte Ausnutzung der Gradienten-Hardware 
moglich. 

Aufgabe der Erfindung ist es, ein Verfahren zur Steuerung einer Pulssequenz bzw. eine Vorrichtung zur Durchfiihrung 
des Verfahrens so auszugestalten, daB die obengenannten Nachteile vermieden werden. 

ErfindungsgemaB wird diese Aufgabe durch ein Verfahren nach Anspruch 1 bzw. eine Vorrichtung nach Anspruch 8 
45 gelost. Hierbei wird aus Sequenzdaten, die im k-Raum vorgegeben sind, wahrend der Laufzeit der Pulssequenz ein Steu- 
erdatensatz fur Gradientenverstarker und einen Hochfrequenz-Sende- und Empfangskanal berechnet. Darnit ist ein Ein- 
griff in die laufende Sequenz, also eine Anderung der Sequenzparameter moglich. Durch eine geringe Latenzzeit von 
maximal einer Zeitscheibe ergibt sich ein geringer Speicherplatzbedarf. 

Vorteilhafte Ausgestaltungen der Erfindung sind in den Unteranspriichen angegeben. Ein Ausfuhrungsbeispiel der Er- 
50 findung wird nachfolgend anhand der Fig. 1 bis 11 naher erlautert. Dabei zeigen: 

Fig. 1 schematisch den Aufbau eines Kernspintomographiegerats, 

Fig. 2 etwas detaillierter den Aufbau der Anlagensteuerung, 

Fig. 3 einen Teil einer typischen Pulssequenz, 

Fig. 4 die Belegung einer Rohdatenmatrix im k-Raum 
55 Fig. 5 ein FluBdiagramm einer herkommlichen Steuerung, 

Fig. 6 die herkommliche Einteilung einer Pulssequenz in einzelne Zeitscheiben, 

Fig. 7 die herkommliche Definition einer Gradientenanstiegsflanke, 

Fig. 8 das genereile Konzept der erfindungsgemaBen Steuerung. 

Fig. 9 die Aufteilung der Sequenzsteuerung, 
60 Fig. 10 ein FluBdiagramm einer Ausfiihrungsbeispiels der neuen Steuerung, 

Fig. 11 ein Gradientendiagramm. 

In Fig. 1 ist stark schemausiert der ubliche Aufbau eines Kernspintomographiegerats dargestellt. Dieses besteht aus ei- 
nem ein homogenes Grundfeld erzeugenden Magnetsyslem 1 bis 4, das von einer Stromversorgung 11 gespeist wird. Im 
Magnetsystem sind Gradientenspulensysteme 7, 8 vorgesehen, die von einem Gradientenverstarker 12 angesteuert wer- 
65 den. Die Gradientenspulensysteme sind zur Erzeugung von Magnetfeldgradienlen in drei Raumrichtungen x, y, z eines 
Koordinatensystems 6 ausgefuhrt. Das Untersuchungsobjekt 5 ist von einer Hochfrequenzantenne 9 umgeben, die mit ei- 
ner Hochfrequenz-Sendeeinheit 14 sowie mit einer Hochfrequenz-Empfangseinheit 15 verbunden ist. Die Hochfre- 
quenz-Sendeeinheit 14 und die Hochfrequenz-Empfangseinheit 15 sind Bestandteil eines Hochfrequcnzsystems 16, in 
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dem unter anderem die empfangenen Signale abgelastel und phasenempfindlich demodulien werden. Aus den demodu- 
herten S.gnalen w,rd mil c.nem Bildrekonstruklionsrechner 17 ein Bild erstellt, an einen Hostrechner 20 weitergegeben 
und auf emem Monitor 18 abgebildet. Die gesamte Einheh wird von einem Steuerrechner20 angesteuen Die Hardwarc- 
Steuerung ubermmmt ein Steuerungsrechner 21, der an den Hostrechner 20 angeschlossen isi. Die gesamte Anlagen- 
steuerung ist in Fig. 2 noch etwas genauer dargestellu Der Steuerungsrechner 21 enthalt eine Standard-CPU auf der die 
eigenthche MeBsequenz ablauft. An diese CPU sind mindestens eine Hochfrequenz-Sendeeinheit 14, mindestens eine 
Hochfrequenz-Empfangseinheit 15 und drei digitale Signalprozessoren (DSP) 22, 23 und 24 gekoppelt. Die Gradienten- 
DSPs konnen via Interrupts Daten von der Steuerungs-CPU anfordem, sobald ihr lokaler Puffer leereelaufen ist Weiter 
besleht eine digitale Verbindung zum Bildrekonstruklionsrechner 17, um die Daienerfassung mil der Bildrekonstruktion 
zu synchronisieren. Ein Bus zum Hostrechner 20 erhalt eine Echtzeit- Verbindung zwischen dem Host und der Steue- 
rungs-LPU aufrecht, um dem Benutzerdie Moglichkeit zu geben, auch wahrend der Messung auf den MeBablauf einzu- 
wirKen. 

Ein Beispiel fur eine Pulssequenz ist in Fig. 3 dargestelli. Jede Pulssequenz lafit sich einteilen in eine Sendephase in 
der - unter Umstanden unter gleichzeitiger Einwirkung eines definierten Gradienten (Gs in Fig. 3) - ein Hochfreouenz- 
puls emgestrahlt wird. SchlieBlich gibt es den "WARP'-Zustand, in dem z. B. durch ein definiertes Gradientenzeitinte- 
gral eine Phasencodierung der Kemspins erfolgt oder auch eine spezieUe FluBcodierung. eine bestimmte Diffusionsco- 
dierung usw. . Der WARP-Zustand dient femer dazu, die Kontinuitat und Stetigkeit der Gradientenstrome zwischen den 
angrenzenden Zeitscheiben herzusteUen. Wahrend des WARP-Zustands sei keine Hochfrequenzaktivitat erlaubt 

Weiter gibt es den Empfangszustand (in Fig. 3 mit ADC bezeichnet), in dem mit einem Analog/Digital-Konverter das 
entstehende Signal abgetastet wird. Dabei isl die Moglichkeit eines gleichzeitig definierten Feldgradienten (GR in Fie 3) 
vorgesehen. b ' ' ' 

SchlieBlich gibt es Wartezustande, in denen keine Gradientenfelder akuv sind und weder gesendet noch empfansen 
wird. Diese konnen als Sonderfall des WARP-Zustands gesehen werden. 

Die Pulssequenz nach Fig. 3 wird N mal mit unterschiedlichen Amplituden-Zeitintegralen des Phasencodiercradien- 
ten GP wiederholt, so daB man N unterschiedlich phasencodierte Signale S erhalt. Jedes Signal wird m mal abgetastet 
und die.emzelnen MeBwerte in eine Zeile einer in Fig. 4 schematise dargesteUten Rohdatenmatrix RD eingetragen 
Diese Rohdatenmatnx kann man als MeBdatenraum betrachten, der im allgemeinen als k-Raum bezeichnet wird Fiir den 
k-Kaum gilt folgende Definition: 
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k x (t) = yj G x (f )dt" 
0 



k y (t) = yj G y (f )dt' 

0 

t 

k z (t) = yj G 2 (t' )df 
o 

Dabei ist y die Larmorkonstante und G„ G,, G 2 ein Magnetfeldgradient in der Richtung x, y bzw. z eines kartesischen 
Koordinatensystems. Auf die Pulssequenz nach Fig. 3 ubertragen, konnte z. B. der Schichtselektionsgradient GS in z- 
Richtung der Phasencodiergradient GP in y-Richtung und der Auslesegradient GR in x-Richtung Uegen, so daB fiir die- 
sen >all also gilt: 

G £ = GS,G y =GP,G» = GR. 

Aus einem Rohdatensatz im k-Raum, also der Rohdatenmatrix RD nach Fig. 4, kann man nun ein Bild rekonstruieren 
da zwischen dem Ortsraum (also dem Bild) und dem k-Raum mathematisch der Zusammenhang iiber folgende mehrdi- 
mensionale Fourier-Transformation besteht: 

f rr . i(k x x + k v y + k_z) 

S(k x , k y/ k z ) = //J p(x, y, z)e x ■ y* z dxdydz 

Dabei ist p(x,y,z) die Spindichteverteilung und S das erhaltene Signal. Da die MeBwert e als diskrete numerische Werte 
vorliegen wird die Founer-Transformation als diskrete Fourier-Transformation mittels FFT (Fast Fourier Transform)- 
verrahren durchgefiihn. 

Wie bereiis oben erwahnt, belegt jedes Signal S eine Zeile der Rohdatenmatrix. Die Zeilenposition ist dabei entspre- 
chend den ob.gen Ausfuhrungen durch den Wert ky, also durch das Zeitintegral iiber den vorausgehenden Phasencodier- 
gradienten GP fesigelegi. 

Aufgabe der PulssequenzsLeuerung isi es nun, Hochfrequenzsender, EmpfangskanaJ und Gradienlenverstarker so zu 
steuern, daS genau die gewunschie Pulssequenz erzielt wird. Dabei gibt es eine Reihe von unterschiedlichen Sequenzty- 
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pen, die man alien als Sequenztopologie bezeichnen kann. In einer grobcn Einteilung sind dies z. B. Gradicntenechose- 
quenzcn und Spinechoscquenzen. Dicse Gruppen von Pulssequenzen lassen sich wiederum feiner unterteilen, z.B. in 
FLASH -Sequenz^n. FlSP-Seqiienzen. Turbospinechosequenzen und EPI-Sequenzen, urn nur einige Beispiele zu nen- 
nen. Urn die Bedienung so einfach wis ntogiieh zu n:ach£n r slellt man in einem Menu vcrschiedene anatomische Berei- 

5 che (z. B. lumbale Wirbelsaule) mil jeweils optimierten Pulssequenzen zur Verfugung. Der Benutzer sol) jedoch anderer- 
seits die Mogkchkeit haben, die Pulssequenzen auf seine speziellen Bedurfnisse zuzuschneiden. 

Ein FluBdiagramm fiir den Ablauf einer herkommlichen Sequenzsteuerung ist in Fig. 5 dargestellt. Der Benutzer kann 
zunachst ein MeBprotokoll (z. B. lumbale Wirbelsaule) aus einem Menu auswahlen. Ein HOST-Computer selektiert dar- 
aufhin einen Standard fur eine Pulssequenz, die fur diesen anatomischen Bereich in der Regel die opUmalen Ergebmsse 

10 liefert, und zeigt diese an. Der Benutzer hat aber noch die Moglichkeit, diese Pulssequenz nach semen eigenen Bedurf- 
nissen zu modifizieren. Zu jeder Pulssequenz gehoren noch verschiedene Parameter, z. B. die RepetiUonszeit, die Echo- 
zeit das Field of View, die Schichtdicke, die MatrixgroBe und der Flipwinkel. Der HOST-Computer wahlt aufgrund der 
Anatomie-Selektion auch fur diese Parameter einen Standard aus und zeigt diesen an. Der Benutzer hat auch hier die 
Moglichkeit selektierte Parameter zu modifizieren und nicht standardmaBig vorgegebene Parameter, z. B. den Schicht- 

15 winkel, zu selektieren. . 

Aufgrund des ausgewahlten MeBprotokolls sowie aufgrund der vom Benutzer eingesteUten Parameter und ^equenz- 
Modifikationen wird nun aus einer "Bibliothek" eine Sequenzdatei geladen und Berechnungen in bezug auf Sequenzpa- 
rameter (z B. Schichtdicke) werden durchgefuhrt. Das Sequenzprogramm setzt die MeBparameter in eine Folge von ' lo- 
gischen" Gradientenpulsen urn. Aus der Uberlagerung dieser Pulse entstehen drei lbgische Gradientenstrom-Zeitfunktio- 
20 nen im logischen Koordinatensystem PRS (Phasencodier-, Readout- und Schichtrichtung) Diese werden dann zeitaqm- 
distant diskretisiert in einem Gradientenverlaufsspeicher zusamrnengestellt. Urn Schichten realisieren zu konnen, die m 
einem vom Gradientenspulensystem vorgegebenen physikalischen bzw. kartesischen Koordinatensystem (x,y,z) schrag 
liegen werden die im logischen Koordinatensystem gegebenen Stromsollwerte punktweise in das physikalische Koordi- 
natensystem abgebildet. Hierzu kann man z. B. eine Matrix- Vektor-Multipliziereinheit einsetzen, wie sie im US-Patent 
25 5,349,296 beschrieben ist. . . ■ 

Aufgrund dieser Berechnungen halt man einen Datenstrom von Stromsollwerten, die punktweise an Digital-Analog- 
Konverter zur Ansteuerung der Gradientenverstarker weitergegeben werden. ; 

Wie in Fig 5 gezeigt, wird fur die Ansteuerung die gesamte Pulssequenz in einzelne Zeitscheiben zerlegt, in denen je- 
weils eine Aktion ausgelost wird. Diese Aktion kann z. B. ein Konstanthalten aller Gradienten, eine Anstiegsflanke ernes 
bestimmten Gradienten und eine Abfallflanke eines bestimmten Gradienten sein. Die Gradientenflanken selbst werden 
definiert durch einen Startwert und eine Anderung A innerhalb der Zeitscheibe. Der Verlauf der Gradientenflanke ist in 
einer Tabelle abgelegt, auf die die Pulssequenzsteuerung zugreif t. Diese Tabelle kann - wie in Fig. 6 dargestellt - erne ii- 
neare Gradientenanstiegsflanke definieren. Es konnen aber auch noch andere Tabellen fur Gradientenflanken, z. B. eine 
sinusforrnige Flanke, vorliegen. TT L£ ♦ ^ 

Wie bereits eingangs erwahnt, ist dieses Verfahren, bei dem der Verlauf der Gradienten, Hochfrequenzaktivitaten und 
der Empfangsaktivitaten im wesentlichen vor dem Start der Messung berechnet werden, wemg flexibel, da es keine An- 
derung der Parameter in der laufenden Sequenz ermoglicht. Ferner weist es, da alle vorberechneten Werte abgespeichert 
werden rnussen, einen hohen Speicherplatzbedarf auf. SchlieBlich wird aus den eingangs erlauterten Grunden die mog- 
liche Performance der Gradientenverstarker wegen der notwendigen Worst-Case-Dimensiomerung nicht optimal ausge- 
40 nutzt Die punktweise Abbildung vom logischen in das physikalische Koordinatensystem ist auBerst rechemntensiv. 

GemaB der Erfindung geschieht die Rotation der logischen in die physikalisch realisierten Gradienten-Zeitfunktionen 
durch Abbildung der Gradienten-Zeit-Integrale (GT-Momente) vom logischen in das physikalische Koordinatensystem 
vor dem Entwurf der Gradientenstrom-Zeitfunktionen. Die Gradientenstrom-Zeitfunktionen werden danach zeitschei- 
benweise in Echtzeit aus den Amplituden der angrenzenden Zeitscheiben und den zu realisierenden Gradienten-Zeit-In- 
tegralen fur jede phvsikalische Gradientenachse berechnet und in Form von Gradientenamplituden-Zeitpunkt-Tupeln an 
die digitalen Signalprozessoren fur die Gradienten iibergeben. Dadurch wird die rechenintensive und zeitkrinsche Abbil- 
dung der Stutzstellen-Datenstrome vermieden. Die Drehung der Gradientenzeit-Integrale vom logischen in das physika- 
lische Koordinatensystem sowie die Berechnung der Gradientenamplituden-Zeitfunktion geschieht, in der Steuerungs- 

C ™ Verlauf ei'ner Messung erhalt jeder der drei digitalen Signalprozessoren 22 bis 24 nach Fig. 2 von der Steuerungs- 
CPU 21 eine Kette von Gradientenamplituden-Zeitpunkt-Tupeln, die Eckpunkte im Gradienten-Zeitverlauf darstellen. 
Die digitalen Signalprozessoren 22 bis 24 interpolieren in Echtzeit aufeinanderfolgende Tupel zu Gradientensollwerten 
auf einem diskreten Gradientenabtastraster. Erreicht dabei die Systemuhr die Sollzeit des Zieltupels, so wird dieses ver- 
worfen und durch das nachste Tupel ersetzt, bis auch dessen Sollzeit erreicht ist usw. . Die jeweiligen Zielzeiten rnussen 
55 nicht auf dem Abtastraster liegen, so daB es mogbch ist, die Gradientenrampen mit einer wesentlich feineren als der 
durch das Abtastraster vorgegebenen zeitlichen Auflosung zu verschieben. Damit kann das realisierte Gradientenzeitin- 
tegral auch bei grober Diskretisierung der Gradientenamplitude prazise kontrolliert werden. 

Neben dieser Interpolation errechnen die digitalen Signalprozessoren 22 bis 24 eine Vorverzerning der Gradienten- 
strbme zur Kompensation der gradienteninduzierten Wirbelstrome und Frequenzoffset-Werte zur Kompensation gra- 
60 dienteninduzierter Verschiebungen des Grundmagnetfeldes mittels Verstimmung eines Hochfrequenzsynthesizers, der 
das Hochfrequenzsystern steuert, . 

Das generelle Konzept eines solchen Vorgehens ist in Fig. 8 dargestellt. MR-Sequenzen lassen sich in menrere Hier- 

archiestufen aufgliedem: 

65 - Die auBere Schicht reprasentiert die Sequenztopologie, z. B. FLASH, Spinecho, Turbospinecho usw. 

- Die darunterlieoende Schicht (grobe k-Raum-Struktur) reprasentiert z. B. die Segmentierungsart des k-Raums. 
Beispieisweise gibt es Sequenzen, bei denen der k-Raum nicht in zeitlich linearer Reihenfolge belegt wird oder in 
denen nur der halbe k-Rauin erfaBt wird. 
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- Die innere Schicht (feine k-Raum-Struktur) cnLspricht deiaillierten Problemen, wie der Verteilung eines Phasen- 
codiergradienten auf mehrere Zeitscheiben, z. B. vor und hinter einem 180°- Puis bei Spinechosequenzen. 

- Die innerste Schicht (Einzelpuls-Gradientensynthese) reprasentiert z. B. die Synthese der Gradien ten form zwi- 
schcn zwei Hochfrequenzakiivitaien. 

5 

Nur die innerste Schicht ist hardwarenah, d. h., sie benotigt eine detaillierte Kenntnis der Gerateparameter. Alle auBe- 
ren Schichlen beeinflussen direkt das Konlrastverhalten und reprasentieren nicht maschinenspezifische Kenntnisse, bei- 
spielsweise Relaxations zei ten und Entscheidungen betreffend die Sequenztopologie. 

Die innerste Schicht (und nur diese) stellt eine Absiraktionshulle um die Scanner-Hardware dar und bietet sich daher 
fur eine Automatisierung an. Dementsprechend wird bei der Sequenzsteuerung wie folgt vorgegangen: io 

Die Berechnung des Steuerdatensatzes wird wahrend der Laufzeit der Sequenz vorgenornrnen. Dadurch kann der Se- 
quenzablauf mit rainimaler zeitlicher Latenz nicht deterministisch gesteuert werden. AuBerdern entfallt der Speicher- 
platzbedarf fur den Steuerdatensatz, da nur der gerade benotigte Teil "on demand" berechnet wird. 

Die Berechnung des Steuerdatensatzes wird in zwei Teile geteilt, namlich die oben erlauterte innerste Schicht und den 
Rest. Der der innersten Schicht entsprechende Teil wird vom Betriebssystem der Steuerung geleistet, die anderen Hier- 15 
archiestufen entsprechen dagegen dem eigentlichen Sequenzprogramrn. 

Dieses Vorgehen ist schematisch anhand eines Blockschaltbilds nach Fig. 9 dargestellt Es wird zunachst vom Bedie- 
ner eine Sequenztopologie ausgewahlt und dazugehorige Parameter, wie z. B. Repetitionszeit TR, Echozeit TE, Betrach- 
tungsfenster FOV, Schichtdicke, MatrixgroBe, Flipwinkel und Schichtwinkel. AnschlieBend wird das Sequenzprogramrn 
erstellt. Dieses gibt Spezifikationen in Form von k-Raumvektoren, von Zeitangaben und vom Sende- bzw. Empfangszu- 20 
stand, an die sogenannte k-Raummaschine weiter, die aus diesen GroBen wahrend des Sequenzablaufs Steuersignale an 
Gradientenverstarker und Hochfrequenzsende- und Empfangseinheit befert. Der BegrifF "k-Raum-Maschine" wurde 
deshalb gewahlt, weil hier nicht wie iiblich Gradientenamplituden, sondern direkt k-Raumwerte iibergeben werden, was 
die Programmierung wesentlich vereinfacht. In einem Zustandsmodell kann man die Zustande der Steuerung wie folgt 
beschreiben: 25 

Bei der Kernspintomographie hat man zwei Mittel zur Verfugung, um das Kemspinkollektiv zu beeinflussen: Hoch- 
frequenzpulse und Gradientenfelder. Die Steuerung dieser Mittel erfordert im Prinzip vier Zustande: 

- Wartezu stand 

- Hochfrequenz-Sendeaktivitat 30 

- Empfangsaktivitat 

- Warp- Zu stand. 

Das die Anlagensteuerung treibende Sequenzprogramrn liefert wahrend des Ablaufs der Messung eine Folge von Spe- 
zifikationen von Zeitscheiben, die jeweils einem der drei Zustande angehoren. Die Anlagensteuerung berechnet aus die- 35 
sen Angaben die in der nachsten Zeitscheibe zu realisierende Gradientenstromform und gegebenenfalls Hochfrequenz- 
Aktivitat. Dadurch ergibt sich eine Latenz von maximal einer Zeitscheibe. Die Spezifikation der Zeitscheiben erfolgt ab- 
strakt, z. B. durch direkte Angabe von k-Raumpositionen. Auf diese einzelnen Zustande wurde oben bereits eingegan- 
gen. Besonderer Erlauterung bedarf noch der Warp-Zustand: 

Die optimale Gradientenstromform wahrend des Warp-Zustands kann z. B. vom Steuerungsrechner aus der vom Se- 40 
quenzprogramm gelieferten Angabe der zuruckzulegenden k-Raum-Distanz vollautomatisch berechnet werden. Die vom 
Sequenzprogramrn definierten k-Raurn-Vektoren reprasentieren auch die (u. U. schragen) Schichtpositionen. Es wird zu- 
nachst eine Vektorrotation in das physikalische Koordinatensystem angewandt und erst dann die optimale Gradienten- 
stromform berechnet. Damit vermeidet man das eingangs geschilderte Problem herkornmlicher Anlagen, dafi bei der 
Umsetzung vordefinierter logischer Gradienten in physikalische Gradienten komplexe Stromzeit-Verlaufe entstehen, die 45 
vom Gradientenverstarker nur schwer zu beherrschen sind. AuBerdern kann das Gradientensystem auch bei schragen 
Schichten besser ausgenutzt werden, da man nicht - wie bei herkornmlichen Anlagen - bei der Sequenzdefinition vom 
worst case beziiglich der Slewrate und der Amplitude ausgehen muB. SchiieBlich wird bei diesem \forgehen die Program- 
mierung der Sequenz erheblich vereinfacht. Sowohl die obengenannte Vektorrotation als auch die Berechnung der opti- 
malen Gradientenstromform erfolgt wahrend des Sequenzablaufs. 50 

In Fig. 10 ist ein Blockdiagramm fur den Zeitablauf der Steuerung vorgesehen. Wie beim herkornmlichen Verfahren 
nach Fig. 5 wahlt der Benutzer zunachst ein MeBprotokoll aus und der HOST-Computer selektiert dann verschiedene 
Standardwerte, die der Benutzer modifizieren kann. An diesem Punkt startet dann die Messung, wobei wahrend des Se- 
quenzablaufs der Steuercomputer aufgrund der Vorgaben Sollwerte fur Hochfrequenz-Sendeeinheit und Gradientenver- 
starker berechnet und den ADAC-Wandler steuert. 55 

Anhand von Fig. 11 wird im folgenden die Berechnung eines optimalen trapezformigen Dephasier-Gradientenpulszu- 
ges mit einem bestimmten Strom-Zeit-Integral dargestellt 

Bei konstanter Slewrate G und minimaler Anzahl Gradientenrampen ergibt sich fur einen Dephasierpuls nullter Ord- 
nung bei einem Anfangsgradient G und einem Endgradient G 2 bei einer Gesamtzeit T und einem Gradienten-Zei [integral 
(k-Raum-Distanz ° 60 

T 

M 0 =J G dt 
0 

65 

der unbekannte Zwischengradieni G\. 
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M 0 - i (G 0 + G 2 )t 



G ._ i ..- (i; 

1 T - t 



G 2 -Gq] 

mit t = J : . (2, 

Falls damit Gi > max(G 0 , G 2 ), ist Gi tatsachiich 
G { = max(G 0 ,G 2 ) + rG (3). 

15 

Entsprechend fur Gi < min(Go, G 2 ), 
G { = min (G 0 ,G 2 ) - rG (4) 
20 mit dem Hilfsausdruck 



T - t IT 2 - t 2 



r = - 



M 0 " g (G ° + G2)T 



25 2 



(5) 



Falls das Wurzelargument negativ ist oder Gi das Amphtudenlimit uberschreitet, ist der gesamte Gradientenpulszug 
innerhalb der Spezifikationen nicht zu realisieren. 

Die Rampenzeiten to! und tj und die Dachzeit tj ergeben sich aus 



30 



35 



tl2 = El^l (7) 
G 

40 t^T-to!-^ (8) 

Patentanspriiche 

1. Verfahren zur Steuerung der Pulssequenz fur eine MR-Anlage, dadurch gekennzeichnet, daB ein Sequenzpro 
45 grarnm in Form einer k-Raum-Struktur vorgegeben wird und daB ein Steuerdatensatz im physikalischen Bereich fur 

Gradienten, Hochfrequenzpulse und Abtastung von Kernresonanzsignalen wahrend der Laufzeit der Pulssequenz 
berechnet wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB ein die Steuerung treibendes Sequenzprogramm wah- 
rend der Laufzeit der Pulssequenz eine Folge von Spezifikationen von Zeitscheiben liefert, die jeweils einem der 

50 drei Zustande 

- Wartezustand 

- Hf-Akuvitat (Senden oder Empfangen), ggf. gleichzeitig rnit einem definierten Gradienten 

- Warp-Zustand 

angehbren und daB aus diesen Angaben wahrend der Laufzeit der Sequenz die in der nachsten Zeitscheibe zu 
55 realisierende Gradientenstromform und ggf. Hf-Aktivitat errnittelt wird. 

3. Verfahren nach Anspruch 1 oder 2, wobei Gradienten-Zeit-Integrale aus einem logischen in ein physikalisches 
Koordinatensystem transformiert werden und wobei dann Gradientenstrom-Zeitfunkuonen gebildet werden. 

4. Verfahren nach Anspruch 3, wobei Gradienten-Zeitfunktionen zeitscheiben weise in Echtzeit aus den Gradien- 
ten amplituden angrenzender Zeitscheiben und den in der betreffenden Zeitscheibe zu realisierenden Gradienten- 

60 Zeit-Integral berechnet werden. 

5. Verfahren nach Anspruch 4, wobei Gradientenstromsollwerte durch eine Folge von Gradientenamplituden-Zeit- 
punkt-Tupeln vorgegeben werden, die Eckpunkte im Gradienten-Zeitverlauf darstellen, wobei in Echtzeit aufeinan- 
der folgende Tupel zu Gradientenstromsollwerten auf einem diskreten Gradienten-Abtastraster interpoliert werden. 

6. Verfahren nach Anspruch 5, wobei ein Ziel tupel der Interpolation verworfen und durch das nachste Tupel ersetzt 
65 wird, sobald eine Systemuhr der Steuerung die Sollzeil des Zieltupels erreicht hat. 

7. Verfahren nach einem der Anspriiche 1 bis 6, dadurch gekennzeichnet, daB die Spezifikation des Warp-Zustands 
fur die Steuerung durch Vorgabe von k-Raunivektoren erfolgt. 

8. Steuerung fur eine MR-Anlage mit einem Steuerrechner, der aus einem in Form einer k-Rauiustruktur vorgege- 
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benen Puissequenzprogramm wiihrend der Laufzeit der Pulssequenz einen Steuerdatensatz fur Gradientenversiarker 
und einen Hochfrequenz-Sende- und Empfangskanal berechnet. 
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